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ABSTRACT 

We present results from a multi-year monitoring campaign of the broad-line 
radio galaxy 3C 120, using the Rossi X-ray Timing Explorer (RXTE) for nearly 
five years of observations. Additionally, we present coincident optical monitoring 
using data from several ground-based observatories. Both the X-ray and optical 
emission are highly variable and appear to be strongly correlated, with the X-ray 
emission leading the optical by 28 days. The X-ray power density spectrum is 
best fit by a broken power law, with a low-frequency slope of —1.2, breaking to 
a high-frequency slope of —2.1, and a break frequency of logz/b = —5.75 Hz, or 
6.5 days. This value agrees well with the value expected based on 3C 120's mass 
and accretion rate. We find no evidence for a second break in the power spectrum. 
Combined with a moderately soft X-ray spectrum (F = 1.8) and a moderately 
high accretion rate (m/rhEdd ~ 0.3), this indicates that 3C 120 fits in well with 
the high/soft variability state found in most other AGNs. Previous studies have 
shown that the spectrum has a strong Fe Ka line, which may be relativistically 
broadened. The presence of this line, combined with a power spectrum similar 
to that seen in Seyfert galaxies, suggests that the majority of the X-ray emission 
in this object arises in or near the disk, and not in the jet. 
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Introduction 



Variability has long been recognized as a universal property of active galactic nuclei 
(AGN). Produced close to the nucleus, X-ray emission is often the most variable of any 
waveband. Early analysis of X-ray variability showed that the power density spectrum 
(PDS) follows a power law, or "red noise" form, P (x u", where P is the variability power at 
frequ ency u, and a is the power law slope (IMcHardy fc Czernylll987l : iLawrence fc Papadakis 
19931 ). No break frequency was seen, due to the short nature of the observations. Later results 
demonstrated that most AGN can be modeled by a broken power law PDS, where the slope 
changes from beir ig relatively shallow at low frequencies to some steeper value above a critical 
break frequency Jfidelson fc Nandral[l999l : luttlev et allbood : iMarkowitz et allboosh . 



This behavior closely mirrors what is seen in galactic black hole X-ray binaries (BHXRBs), 
which contain black holes with masses of only a few Mq. Much of the work connecting AGN 
and BHXRBs has focused on Cygnus X-1, which is one of the brightest and most well-studied 
BHXRBs. Cygnus X-1 is commonly in one of two different variability states: a low/hard 
state with a low flux level, hard X-ray spectrum, and a PDS with 2 separate breaks; and 
a hig h/soft state, with higher X-ray flux, soft (stee p) X-ray spectrum, and a singly-broken 
PDS (IRevnivtsev et al.ll200d : iMcHardy et al.ll2004l ). Based on scaling of the break frequen- 
cies with black hole mass and accretion rate, nearly all of the AGN studied to date match 
best with the high/sof t state, showing only one break in their power spectrum (Akn 564 is 
the sole exception, see IMcHardy et al.ll2007l ) 



The break frequency scales over more t han 6 orders of magnitu de in mass, if one extrap- 
olates the break frequencies of Cygnus X-1 (IMarkowitz et al.ll2003l ). The variability state is 
thought to be driven by the accretion rate, with the high/soft state correl ating to a higher 
accretion rate which may push in the inner edge of the accretion disk (lEsin et al.l 119971 : 
Churazov et al.| 2001 ). This has been combined into a 'fundamental plane' of variability by 
McHardy et al.l (120061 ) , with mass and bolometric luminosity uniquely determining the break 
timescale. 

In Seyfert 1 galaxies, where no jet is detected, the X-ra.y eni ission is mainly thought to 



arise in or near the accretion disk (iHaardt fc Maraschi 



1991 



1993h. T his variability may also 



be connected to the broad line region (BLR), as IMcHardy et al.l (120061 ) showed that the break 
timescale is also related to the full-width half-maximum (FWHM) of H/3 emission. Thus far, 
efforts to analyze the X-ray variability of AGN with jets have been limited. Furthermore, 
the connection between the X-ray variability and the jet in AGN needs to be defined better 
than has been possible thus far. 



The nearby {z = 0.033) broad-line radio galaxy (BLRG) 3C 120 is an interesting target 
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for this type of variability analysis. Originally clas s ified as a type 1 Seyfert, 3C 120 was 



first discovered as an X-ray source by iForman et al.l (11978 ) using the Uhuru observatory. A 



superluminal radio jet was later disc overed by 



extreme optical variability at times (IWebb et al 



Walker et al 



1988 



fll982h. 3C 120 h as also exhibited 
1980[). with variations 



French fc Miller 



of more than half a magnitude on timescales of a few days. iMarscher et al.l (120021 ) used 
coordinated VLBI and RXTE observations to show that dips in the X-ray flux corresponded 
with ejections of superluminal knots of bright emission along the jet. 

The structure of our paper is as follows: in §2 we discuss our observations and data 
reduction methods, discuss the stationarity of the light curve in §3, present results for X-ray 
and optical correlation in §4, present results from Monte Carlo simulations of the power 
spectrum in §5, and in §6 we discuss implications of our results and present conclusions. 



2. Observations and Data Reduction 

We have observed 3C 120 for several years at both optical and X-ray wavelengths, on a 
variety of timescales. In Table 1 we summarize our monitoring at each timescale. Below we 
detail our data reduction methods. 



2.1. X-ray Data 

All of our X-ray observations come from a multi-year monitoring program using RXTE. 
3C 120 was monitored three times per week from 2002 March until 2007 January, excluding 
several 8-week periods when the object's proximity to the Sun prevented observations. We 
refer to these data as the 'long term' light curve. 

Additionally, from 14 November 2006 until 12 January 2007 3C 120 was observed every 
6 hours during this period. These data are referred to as the 'intermediate' sampling data. 
Finally, 3C 120 was monitored nearly continuously from 13-22 December 2002, which we 
refer to as the 'short' term light curve. The short term light curve was binned to an interval 
of 4000s, because AGN typically do not show variability on shorter timescales, and such 
binning also improves computational timescales greatly. 

We use only data taken by PCU detectors and 2, in STANDARD2 data mode. All 
of our data were reduced using FTOOLS v5.2 software, provided by HEASARC. Data were 
excluded if the Earth elevation angle was < 15°, pointing offset > 0.02°, time since South 
Atlantic Anomaly passage < 30 minutes, or electron noise > 0.1 units. Counts were extracted 
from the top PCU layer only to maximize the signal to noise ratio. 



-4- 



Figure 1 shows the hght curves for all three timescales. In all subsequent analysis, we 
do not interpolate any gaps in the light curves with the exception of the short term light 
curve, which has a larger number of gaps due to Earth occultations and other observing 
programs, and a lower variance than the other data sets. The 2-10 keV X-ray spectral 
index remained relatively constant throughout our observations, with all ObsIDs showing an 
average Fx = 1.8 with a standard deviation of 0.1, where vF^ cx 



2.2. Optical Data 

We have also monitored 3C 120 at optical wavelengths during most of our long term 
X-ray observations using multiple ground-based observatories. From 2004 September until 
2007 January, optical data were taken using the 2 m Liverpool Telescope at La Palma, 
Canary Islands, Spain. We also supplement these observations using data from the 1.8 m 
Perkins Reflector at Lowell Observatory in Arizona, and the 1.3 m SMARTS consortium 
telescope at Cerro Tololo Inter-American Observatory in Chile. Data were reduced using 
standard methods, and a correction applied, if necessary, to match the flux measured by the 
Liverpool Telescope. 

All of our optical data were binned to the same 2.33 day timescale as our long term 
X-ray light curve, giving us a total of 124 data points. Figure 2 shows the binned optical 
light curve, along with the relevant portion of the X-ray light curve. 



3. Stationarity 

Before we can undertake any meaningful analysis of variability, we must first determine 
whether or not the data are stationary. In a statistical sense, stationarity means the underly- 
in g variability process doe s not c hange over time. We can test for this using the "5"' statistic 



of iPapadakis fc Lawrence! (119951 ). which measures stationarity in the power spectrum. The 
light curve is divided into two equal halves, the power spectrum is computed for each half, 
and then the two are compared using the S statistic, where S is defined as: 

g^^^ l0gPi(z/)-l0gP2M 



A/var[logPi(z/)] + var[logP2('^)] 



where P|,2(^) are the power s pectra for each half of the light curve, and var[logPi_2] = 0.31 
(IPapadakis fc Lawrencelll993l ). We can then sum over all frequencies to find the test statistic 
5, 
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n=l 

where N is the total number of frequencies for which we have calculated the power spectrum. 
For two identical power spectra, S will have a Gaussian distribution with zero mean and 
unit variance. Or in other words, if 1 5*1 < 1 we can reasonably assume the variability is 
stationary. 

For our long-term light curve, we find that 5* = 0.238. This indicates that the variability 
in 3C 120 is stationary, and that the fundamental variability mechanism does not change 
during the 4+ years of our observations. 



4. X-ray /Optical Cross Correlation 



To examine the relationship between the X-ray and optical emission, we use the cross 
correlation function (CCF). The traditional CCF requires evenly sampled data, and can be 
computationally intensive. Because of observing constraints, neither our X-ray nor optical 



d ata are evenly samp 



ed. T o solve this issue, we use the discrete correlation function (DCF) 



of lEdelson &: KrolikI (jl988l ). which allows for cross correlation of two unevenly sampled data 
sets. 

The DCF is shown in Figure 3, with positive lag indicating X-ray emission leading the 
optical. Because the DCF has a broad shape, determining a precise maximum value may not 
be the best way to quantify any possible lag, so we fit the DCF with a Gaussian curve. The 
curve provides an excellent fit to the DCF, and has a maximum at a peak lag of -|-28 days, 
at a correlation coefficient of 0.8. 



As discussed in lUttley et al.l (120031 ) . traditional error bars are inadequate for assess- 
ing the significance of the DCF, because adjacent data points in the light curve are "red 
noise" data and not uncorrected. Therefore we use Monte Carlo simulations to assess the 
significance of the correlation. Our method is as follows. 

We began by simula ting 2 independent red noise hght curves, using the method of 
Timmer fc Koenig] (119951 ). One light curve was given a power spectrum identical to the 
X-ray PDS in §5, and the other light curve was simulated using the measured optical PDS, 
which shows a slope of —1.6 with no measured break (Ryle et al. 2009, in preparation). The 
two light curves were then re-sampled in the same fashion as the original light curves, and 
random noise was added in the form of a Gaussian random with mean of zero and standard 
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deviation equal to the average observed error. The hght curves were then correlated, looking 
for any cases where the maximum correlation coefficient was greater than 0.8. Out of 1000 
simulations, this was only the case 2 times. Therefore the correlation coefficient of r = 0.8 
seen in the data is significant at more than 99% confidence. 

To explore the significance of the measu red lag, we use the Flux Randomization / Ran - 
dom Subset Selection (FR/RSS) technique (IWhite fc Peterson! Il994j : iPeterson et al.lll998l ). 
Briefly, random errors are added to each light curve, using the same method as above. For 
each light curve with data points total, the same number of data points are chosen at 
random from each curve, and duplicates are rejected. This typically results in a fraction of 
~ 1/e points rejected from each new light curve. The curves are then cross correlated as 
before, the DCF is ffi with a Gaussian, and the centroid noted. This is done a large number 
of times (1000 simulations), and the distribution of DCF centroid values is used to estimate 
uncertainties. 



As noted by iMaoz fc Netzerl (119891 ). distributions of lags are nonnormal and the stan- 
dard deviation is not a good estimate of the uncertainty in the lag. Therefore similar to 
Peterson et al.l (119981 ) . we quote uncertainties at iArgs, where 68.27% of the realizations 
yield results between Tmedian — Atqs and Tmedian + Args, which corresponds to a Icr error for 
a normal distribution. Using this method, we estimate the lag for 3C 120 as 28.73I5 §7 days. 
The distribution of lags measured using the FR/RSS method is shown in Figure 4. 



5. X-ray Power Density Spectrum 



The large amount of X-ray data, covering a wide variety of timescales, allow us to calcu- 
late the power density spectrum (PDS) over more than 4 decades of frequency. Calculating 
the PDS for any discretely sampled, finite length data set is inherently problematic, suffering 
from a number of complications such as windo wing, aliasing , and r ed noise leakage. To avoid 
these complications, we follow the method of lUttley et al.l (120021 ). which we outline briefly 
below. 



We begin by calculating the square modulus of the discrete Fourier transform for each 
light curve, using the formula 



N 



-I 2 



^f{ti) cos(27rz/ti) 



4 = 1 



+ 



TV 



-T 2 



^f{ti) sin(27rz/ti 



i=l 



where ti are the light curve data points, and |F(z/)p is calculated at evenly spaced frequencies 
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from z/min = 1/T to z/Nyq = N/2T, where T is the total duration of t he hght curve. The PD S 



is then normahzed using the fractional rms squared normalization (iMiyamoto et al.lll99ll ). 

where is the average count rate of the light curve. The observed PDS is then binned 
geometrically every factor of 1.5 (0.18 in the logarithm). 

Monte Carlo simulations are then used to find the best-fit model for the observed PDS. 
A simulate d light curve is gen e rated based on an input model power spectrum, based on the 



method of iTimmer fc Koenigl (119951 ). The simulated light curve is created with at least 100 
times the duration of the actual light curve, to avoid the problem of variability power from 
longer timescales being shifted to shorter timescales (red noise leakage). The simulated light 
curve is then split into 100 segments of equal length. 

Because our data are not continuously sampled, power above the Nyquist frequency can 
also be aliased back into the power spectrum. To avoid this, our simulated light curves have 
a time resolution of O.lTgamp, where Tgamp is the sampling timescale of the light curve. 

The simulated light curves are then re-sampled in the same fashion as the observed light 
curves. The PDS is then calculated for each of the 100 curves, and binned by factors of 1.5. 
The 100 individual power spectra are then averaged, and error bars for each data point are 
calculated from the rms spread of the individual PDS points. 

Poisson noise is added to the simulated power spectra using 



-Ppoisson 



iV('^Nyq 



where a{i) are the error bars for each point in the light curve. We can also estimate analyt- 
ically the amount of power which will be aliased back into the PDS from frequencies above 
the minimum resolution of the simulations. This power is given by: 

I /-l/CSToxp) 
^Alias = / P{u)du, 



where Tgxp is the exposure time of our observations, typically ~ 1 ksec. 

To determine the goodness of fit of the model, we use a modified fit. First we compare 
the observed PDS to average of the simulated power spectra: 
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Xobs 



(AP„„(i/))= 



where Pobs(i^) is the observed PDS at frequency u, Psim{i^) is the average of the simulated 
power spectra at that same frequency, and APsiin(i/) is the rms spread of the power spectra. 

We then randomly choose 1000 combinations of simulated power spectra at each sam- 
pling rate. We then calculate x^isf where: 



2 V- (PsM - PsMY 



Xd: 



ist 



E 



APsim(i^)' 



Of these 1000 combinations, the fraction where Xobs is smaller than x^ist represents 
the probability of acceptance of the model, i.e., the probability that the model provides an 
acceptable fit to the data. This process is repeated for a grid of slopes and break frequencies, 
to determine the parameters which yield the highest probability of acceptance. 



5.1. Unbroken Power Law 

We began by fitting the observed PDS with an unbroken power law of the form P{i/) oc 
u". This yields a best-fit power law slope of a = —1.5, with a probability of acceptance of 
only 1.3%, and is clearly not an adequate fit to the data. 



5.2. Broken Power Law 

Visual inspection of the PDS indicates that the shape may be a broken power law, with 
a fixed slope at low frequencies that breaks to a steeper slope at high frequencies. For our 
simulations, we use a broken power law of the form: 



where A is the normalization, i\, is the break frequency, ckl is the low-frequency slope, and 
an is the high-frequency slope. The low frequency slope was allowed to vary from —1.0 to 
— 1.5 in increments of 0.1, and the high frequency slope was allowed to vary between —1.5 and 
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—3.5 in increments of 0.2. The break frequency was allowed to vary between logz^b = —8.0 
and —4.0, incrementing the break by factors of 2 (0.30 in the logarithm). Once a fit was 
obtained, a finer grid with increments of 0.1 in slope and 0.0414 in break (factors of 1.1) 
were used to improve the best-fit values. 

Figure 5 shows the observed PDS fitted with the best-fit broken power law model. 
The best fit parameters show a low-frequency slope of = —1.2, high-frequency slope 
ttH = —2.1 ± 0.4, with a break at a frequency of logz/b = —5.75 ± 0.43 Hz, or a timescale 
of 6.51^4^° days. The probability of acceptance of this model is 84.5%. Figure 6 shows 
confidence contours for variations in the break frequency and high frequency slope. Error 
bars for the slope and break were obtained by taking a 1-dimensional slice through the 
confidence contours at the best-fit value, fitting the contours with a Gaussian and finding 
the standard deviation. 



5.3. Doubly-Broken Power Law 

Given that the low/hard state seen in BHXRBs shows two break frequencies, we also 
model the PDS of 3C 120 with a doubly broken power law. For computational efficiency, 
the slopes were fixed at values of 0, —1, and —2, similar to that seen in Cygnus X-1 



( iRevnivtsev et al.ll2000l ). The low frequency break was allowed to vary from logz^b.io = —9.0 
to —5.5, and the high frequency logi^b.hi allowed to vary between 5 and 10,000 times greater 
than the low-frequency break. 

A best fit was obtained at logz/b,io = —8.03, logz/b,hi = —5.83, and a confidence level of 
63.2%. The low frequency break is near the lowest measured frequency for our PDS, effec- 
tively turning the model into a singly-broken power law. Therefore we find no evidence for 
a doubly-broken power law in the PDS of 3C 120, with a singly-broken power law providing 
an adequate fit. 



Discussion and Conclusions 



We have shown that both X-ray and optical emission are highly variable over a period 
of several ye ars, with fractional variabili ties of 21.6% and 8.2% in the X-ray and optical, 
respectively (IMarkowitz &: Edelsonll2004l ). The variations are strongly correlated, with the 
X-ray variations leading the optical by 28 days. Additionally, the X-ray PDS is best fit by a 
broken power law, with a slope of —1.2 breaking to a slope of —2.1 above a break frequency 
of z/b = —5.75 Hz, or 6.5 days. 
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Given that the PDS is best fit by a singly-broken power law, and does not show any 
signs of a double break, 3C 120 appears to fit well with the high/soft state power spectra 
found in most other AGN. This i s appropriate, gi ven the moderately high accretion rate 
of 30% Eddington for this object (lOgle et al.ll2005l ). The X-ray spectrum is also relatively 
steep, with an average photon index of F = 1.8 during our long-term monitoring, similar to 
the values of F = 2 — 3 seen in high/soft state BHXRBs. 



Using the fundamental plane of AGN variability (IMcHardy et al.l l2006l ). we can in- 
vestigate whether the break timescale measured for 3C 120 fits with a relation based on 
radio-quiet objects. The expected break timescale Tb is: 



logTe = 2.1 log Mbh - 0.98 log Lboi - 2.32 

where Mbh is the mass of the black hole in units of lO^ Mf;^, and Lboi is th e bolometric 
luminosity in units of 1 0^^ erg/s. For 3C 120, Lu^^ = 2 1.87 JWoo fc Urrvll2002h . and Mbh = 



SOlls^, as measured by IVestergaard fc PetersonI (120061 ) using reverberation mapping in the 
optical and UV. Using these values, we would expect to find a break at a timescale of 
Tb = 0.86+^:^^ days, with the errors based on the large uncertainty in the mass. Therefore 
our measured break timescale agrees well with the prediction based on the fundamental 
variability plane, within the errors on both the mass and our measured break. 

We can also check the relationship between H/? FWHM and Tb for 3C 120. This 
relationship is given by: 



logTB = 4.2 log H« FWHM - 14.43 



where H/j fwhm is in units of km/s (IMcHardy et al.ll2006l ). For 3C 120. lMarziani et al.l (120031 ) 
measured a FWHM of 2328 km/s. This gives us a predicted break timescale of 0.51 days, 
different by a factor of ~ 2 from the mass/luminosity scaling relationship, and different by 
a factor of ~ 13 from our actual measured value. The HjS — Tb relationship is based on 
a sample of radio-quiet objects. Given that 3C 120 is radio loud, with a superluminal jet, 
the broad-line region may not share the same characteristics as a radio-quiet object. More 
observations of radio-loud AGN are needed to determine if a separate HP — Tb relationship 
exists for these objects. 

Given that the X-ray behavior of this object appears similar to that from radio-quiet 
Seyfert galaxies, we can examine whether the X-ray emiss i on com es from the jet, or from a 
source in/near the accretion disk. As noted by lOgle et al.l (|2005l ) . although 2 X-ray sources 
are coincident with the positions of radio lobes, they contribute < 5% of the overall X-ray 
flux. 
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Closer to th e nucleus, we can use sp ectroscopy to determine how much of the flux orig- 
inates in the jet. iBallantyne et al.l (l2004j ) have analyzed a 127 ks XMM /Newt on observation 



of 3C 120. They flnd a neutral Fe Ka line with an equivalent width of 53 eV. iKataoka et al. 



( 120071 ) flnd evidence for a relativistically broadened line using Suzaku. The strength of the 
narrow Fe Ka line indicates that the majority of the X-ray emission has its origin in an 
unbeamed source, such as the accretion disk or corona, as opposed t o the jet. I f the X -ray 
flux were beamed within a narrow cone, reflection would be unlikely. lOgle et al.l (120051 ) flnd 
a half-opening angle for the continuum of > 50°, indicating the narrow jet is not responsible 
for the X-ray emission. 

If the X-ray emission originates near the accretion disk, then where does the 28 day 
X-ray to optical lag arise from? Low polarization in the optical (~ 0.3%, consistent with 
interstellar dust) indicates that the emission is not associated with the jet. Such a lag is 
also too large to arise from light travel time between different regions in the accretion disk, 
however the communication between X-ray and optical regions may propagate at sub-light 
speeds. A more complete optical light curve may give a more precise value of the delay, to 
conflrm the result and better specify the required propagation speed of the variations. 
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Table 1. Summary of Observations 



Duration 


start" 


End" 






Long 


52335.42 


54114.39 


1779 


2.33^^ 


Medium 


54054.66 


54114.39 


60 


6h 


Short 


52622.39 


52631.48 


9.09 


1000= 


Optical 


53248.73 


54110.54 


862 


2.33^* 



"Modified Julian Date 
''Length of light curve, in days 
Light curve binning 
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Fig. 1. — RXTE 2-20 keV light curves for long (top), medium (middle), and short (bottom) 
timescales. 
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Fig. 2. — Optical R-band (top) and X-ray (bottom) light curves. 
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Fig. 3. — Discrete correlation function for X-ray and optical light curves. Positive lag 
indicates X-ray emission leading the optical. Dashed line shows a best-fit Gaussian curve, 
with centroid of -|-28 days. See text for a discussion of errors. 
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Fig. 4. — Histogram showing cross-correlation lag centroid versus number of simulation 
results yielding that lag for the FR/RSS method. Bin width is 1 day. Median result occurs 
at a lag of 28.73 days, with uncertainties of -|-At68 = 6.19 and —Ares = 5.87 days. 
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Fig. 5. — Observed PDS for long, medium, and short light curves represented by sohd, 
dotted, and dashed hues, respectively. Monte Carlo simulations for each light curve are 
represented by diamonds, squares, and triangles with error bars. Best fit occurs with the 
low frequency slope fixed at —1.2, high frequency slope of —2.1, and a break frequency of 
logi^b = —5.75 Hz, with a probability of acceptance of 84.5%. 
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Fig. 6. — Confidence contours for a broken power law model with low frequency slope of —1.2. 
Contours are at 5, 10, 25, 50, and 75% probability of acceptance. Maximum confidence of 
84.5% occurs at a high-frequency slope of —2.1 and a break frequency of logi^b = —5.75 Hz. 



